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SECTION  I 


INTRODUCTION 

The  primary  function  of  precision  approach  radars  is  the  detection 
and  tracking  of  landing  aircraft  during  final  approach,  oftentimes 
in  the  presence  of  heavy  rainfall  clutter.  Because  the  operation 
of  such  radars  has  been  restricted  to  X-Band  and  because  the  radar 
cross  section  of  rain  is  an  increasing  function  of  frequency,  techniques 
for  discrimination  of  aircraft  from  rain  clutter  are  of  major  concern 
to  the  system  design  process.  Available  discrimination  techniques 
are  (1)  the  use  of  circular  polarization,  (2)  the  use  of  high  range 
resolution,  and  (3)  the  use  of  Doppler  resolution  or  MTI.  For  example, 
the  AN/TPN-19  PAR  uses  a combination  of  all  three,  each  contributing 
toward  the  whole.  However,  this  design  approach  has  proved  costly. 
Currently,  the  AN/GPN-XX  PAR  program  is  considering  approaches  of 
lower  cost  potential. 

In  order  to  provide  an  enlarged  technical  data  base  applicable 
to  the  AN/GPN-XX  program.  Project  7080  entitled  Aircraft  Detection 
in  Rain  was  conducted  in  Fiscal  Year  1974  for  the  Electronic  Systems 
Division,  USAF,  by  the  MITRE  Corporation.  This  project  embraced  three 
specific  tasks:  (1)  the  experimental  measurement  of  the  depolarization 
of  circular  polarization  by  rain  at  X-Band,  (2)  the  investigation  of 
a linearly  polarized,  low  cost,  electronically  scanned  antenna,  and 
(3)  the  application  of  digital  filtering  techniques  to  the  discrimination 


4 


of  rain  clutter  through  Doppler  resolution  of  target  and  clutter. 

As  a result  of  a review  of  this  project  by  the  Rome  Air  Development 
Center  in  its  early  stages,  the  application  of  effort  applied  to 
the  depolarization  measurements  was  emphasized  and  that  applied  to 
the  antenna  and  digital  filtering  tasks  was  diminished. 

This  report  presents  the  results  of  the  digital  filtering 
application  study.  The  results  of  the  study  of  circular  polarization 
degradation  are  given  in  Reference  1 and  those  of  the  low  cost  antenna 
study  are  given  in  Reference  2.  The  principal  conclusions  of  this 
report  are  that  the  suppression  of  rain  clutter  by  40  dB  .solely 
through  the  application  of  digital  filtering  techniques  appears 
feasible  and  that  its  use  in  conjunction  with  the  linearly  polarized 
low  cost  antenna  approach  of  Reference  2 could  possess  significant 
potential  for  a low  cost  system  design.  The  degradation  data  of 
Reference  1 on  the  other  hand  should  be  useful  to  system  design 
concepts  which  utilize  circular  polarization  for  suppression  inasmuch 
as  very  little  data  of  this  kind  is  available. 

In  this  report  digital  filtering  techniques  are  applied  to 
the  specific  goal  of  discriminating  aircraft  approaching  at  represent- 
ative landing  velocities  from  heavy  rainfall  clutter.  This  application 
of  digital  filtering  in  effect  constitutes  a digital  moving  target 
indicator  in  the  sense  that  the  output  is  not  a replica  of  the  input 
signal,  but  rather  a number  related  to  its  amplitude  and  frequency. 

A separate  digitized  number  of  this  type  is  provided  for  each 
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range  resolution  cell.  The  digital  filtering  provides  a band-stop 
characteristic  which  severely  attenuates  the  amplitude  of  fixed  or 
slowly  moving  clutter  and  a band-pass  characteristic  which  does  not 
affect  significantly  the  amplitude  of  targets  moving  at  final  approach 
velocities.  A threshold  can  be  applied  to  the  output  for  automated 
detection  decision  and  range  tracking.  The  technique  utilizes  a 
coherent  un-coded  radar  pulse  of  constant  pulse  duration.  Two  pulse 
repetition  frequencies  are  used  for  blind  speed  elimination.  In-phase 
and  quadrature  detection  are  also  used.  The  physical  implementation 
of  the  digital  circuitry  is  modest  and  is  judged  to  be  within  the 
current  state-of-the-art. 

A further  by  product  of  this  Digital  MTI  is  that  stationary  clutter 
is  suppressed.  It  is  to  be  noted  that  the  target  aircraft  and  the  rain 
clutter  are  in  the  same  range  cell  throughout  this  study. 

Section  II  of  this  report  describes  the  theoretical  basis  for 
Digital  MTI  processing  and  the  application  to  aircraft  detection  in  rain 
clutter.  Section  III  details  the  model  chosen  to  simulate  rain 
characteristics.  Also  described  in  Section  III  is  the  operation  of  the 
filtering  process  upon  sampled  return  signals.  Section  IV  describes  the 
results  of  the  modelling  with  various  filters.  The  results  are  pre- 
sented with  two  conditions.  The  first  is  with  no  aircraft  present, 
the  second  with  aircraft  and  clutter.  Section  V,  with  references  to 
previously  published  results,  delineates  the  limitations  of  the  Digital 
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MTI  due  to  the  rain  clutter's  incoherency  and  to  quantized  noise. 
Conclusions  are  drawn  in  Section  VI  as  to  the  problems  of  coherency, 
and  the  ability  of  filters  to  reject  clutter. 
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SECTION  II 


DIGITAL  MTI  THEORY  AND  DESIGN  APPLICATION 


2.0  SCOPE 

A major  problem  encountered  by  X-Band  precision  approach  radars 
is  the  detection  and  tracking  of  aircraft  targets  in  heavy  rain  fall. 
The  use  of  band-pass  filters  to  separate  the  Doppler  shifted  signal 
of  the  rapidly  moving  target  from  the  clutter  return  of  the  slowly 
moving  rain  offers  a potential  solution.  This  report  presents  the 
results  of  an  investigation  of  the  ability  of  transversal  digital 
filters  to  improve  the  signal- to-rain  clutter  ratio.  The  study  con- 
sisted of  two  parts:  the  development  of  a model  for  the  rain  clutter 
power  density  spectrum,  and  the  implementation  of  a particular  digital 
filter.  Both  parts  were  realized  on  a HP-2100  Digital  Computer. 

The  rain  clutter  was  modelled  statistically  using  realistic 
parameters.  An  input  clutter-to-signal  ratio  of  20  dB  was  modelled. 
This  clutter  corresponded  to  a volume  cell  of  150  m.  x 1°  x 1°  at 

10  Km.  The  radar  cross  section  of  rain  falling  at  50  mm/hr  at  10  GHz 

2 

is  typically  100  m , for  the  rain  cell  discussed  above.  Typical  radar 
cross  sections  of  aircraft  were  considered  to  be  of  the  order  of  1 m^. 

The  digital  MTI  used  was  an  n tap  digital  delay  line  characterized 
by  n pre-determined  weights  w^,  w2»'**wn’  The  tap  spacing  corresponded 
to  the  pulse  repetition  interval.  The  filtering  operation  consisted 
of  multiplying  the  complex  envelope  (amplitude  and  phase)  of  the  com- 
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posite  return  signal  in  a given  range  cell  from  transmission  1 by  w^, 
from  transmission  2 by  W2,  etc.,  adding,  and  then  computing  the  magni- 
tude of  the  sum.  The  results  obtained  with  filters  of  order  n = 15, 

19  and  23  are  presented;  the  subclutter  visibility  improvement  ob- 
tained with  n = 23  is  worthy  of  consideration  for  future  system 
design. 

This  study  was  concerned  with  the  extraction  of  targets  from 
clutter  by  means  of  Digital  Signal  Processing.  The  basis  for  clutter 
rejection  is  that  the  centroid  of  the  Doppler  frequency  distribution 
of  the  clutter  power  spectrum  is  usually  well  below  the  Doppler 
frequencies  of  aircraft  targets.  Thus,  a moving  target  indicator 
(MTI)  Digital  processor  is  envisaged  as  a system  that  has  a stop  band 
characteristic  which  selectively  suppresses  the  clutter  and  a pass 
band  characteristic  which  embraces  the  target  signal. 

The  MTI  processor  considered  in  this  study  could  be  used  with 
a linear  polarized,  step  scanned  radar  of  the  pencil  beam  type. 

This  investigation  concentrated  on  the  processing  of  the  clutter 
and  target  in  one  representative  range-azimuth-elevation  cell.  The 
system  application  of  the  technique  would  employ  parallel  processing 
of  the  other  range  cells,  extending  to  20  Km.  The  aggregate  process 

fel 

would  be  repeated  for  the  different  azimuth  elevation  cells.  L J Of 
necessity,  the  signals  from  each  range-azimuth-elevation  cell  must 
be  processed  in  blocks  of  n pulses. 
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The  composite  radar  return  signals  are  weighted  and  summed.  This 
weighting  operation  corresponded  to  processing  the  signals  through  a 
non-recursive  digital  fitler  and  examining  the  results  of  the  weighting 
process  after  receiving  n pulses.  The  results  of  the  weighting  process 
is  a signal  level  which  is  to  be  compared  with  a pre-determined  thres- 
hold for  detection  of  targets  for  use  as  an  MTI.  This  digital  MTI 
can  also  be  described  equally  as  well  as  a transversal  fitler,  tapped 
digital  delay  line,  n pulse  canceller,  etc. 

2.1  Filter  Characteristics 

The  MTI  processor  design  had  to  satisfy  two  basic  requirements 
for  the  PAR  application:  These  were  (1)  a maximum  range  of  approximately 
20  Km  and  (2)  target  radial  velocities  between  40  and  480  knots, 
approximately. 

A dual  pulse  repetition  frequency  (PRF)  of  8.00  KHz  and  5.33  KHz 
approximates  these  criteria.  The  maximum  unambiguous  range  Rq  as 
determined  by  the  8.00  KHz  PRF  is  given  by 

Ro=  2(5®!'  18-75Km-  (1) 

The  blind  speed,  or  foldover  velocity, is  given  by 

Vf  = X (PRF) /2  (2) 

which  at  3.1  centimeters  is  240  knots  for  8.00  KHz  and  160  knots  for 
5.33  KHz. 
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Figure  1 displays  a general  amplitude  versus  frequency  character- 
istic for  a digital  filter  of  the  type  considered.  It  is  to  be  noted 
that  the  response  of  a digital  filter  is  periodic  at  multiples  of  the 
PRF. 

Figure  2 displays  the  desired  amplitude  versus  frequency  response 
of  filters  operating  at  8.00  KHz  and  5.33  KHz  respectively.  The 
purpose  of  the  5.33  KHz  PRF  and  its  associated  filter  is  to  fill  in 
the  blind  region  (100  to  140  knots)  of  the  8.00  KHz  filter. 

2.2  Operation  of  Digital  MTI 

The  MTI  filtering  consists  first  of  processing  n consecutive 
pulses  of  each  range  bin  at  a PRF  of  8.00  KHz  using  the  8.00  KHz 
digital  filter  and  then  repeating  the  process  at  a rate  of  5.33  KHz 
using  the  5.33  KHz  filter.  Secondly,  whenever  either  output  of  the 
two  filters  exceeds  a predetermined  threshold,  then  a target  detection 
decision  can  be  invoked.  The  actual  setting  of  the  threshold  for  the 
decision  made  upon  the  outputs  of  the  filters  will  of  course  be 
dictated  by  requirements  for  specific  false  alarm  and  detection 
probabilities.  Three  cases  were  studied;  viz,  n = 15,  19,  and  23. 

2.3  Digital  Filters 

In  general,  the  output  response  y at  sample  point  n to  an 
arbitrary  linear  Digital  Filter  can  be  described  by  the  following 
equation: 
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y(n)  = - y(n-l) 


+ a2y(n-2)  + a^  y (n-3)+. . . . affi 


am  y(n-m> 


+ x(n) 


+ 1*2  x(n-l)  + x(n-2)  . . . .w 


^ x(n-r)j  (3) 


Equation  3 states  that  the  response  depends,  in  general,  upon  m 
past  output  responses  and  r past  inputs.  Non-recursive  filters  demand 
that  aq  = a2  = ••  = am  = 0,  which  is  the  case  discussed  in  this  report. 

To  determine  the  frequency  response  of  a sampled  sinusoidal 
frequency  as  input,  the  Z transform  is  performed  on  Equation  3,  and 
e^a^  is  substituted  for  Z.  T is  the  sampling  time,  and  <d  is  the 
angular  frequency  of  the  sinusoid. 


Defining 


00 


(4) 


00 


(5) 


then  the  system  response  function  defined  as  H(Z)  is 


r+1  „ ,-(£-1) 


(6) 
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The  filters  used  in  this  report  had  n elements  or  weights,  thus 
n = r+1,  in  which  case  r=n-l  and  equation  6 becomes 

r H(Z)  = ^ tf  Z’(A-1)  (7) 

£=1 

If  e^"  is  substituted  for  Z in  equation  (7)  then  the  frequency 
response  of  the  filter  is  obtained  for  n samples  taken  at  intervals 
of  time  T.  Table  2 lists  the  actual  frequency  response  of  the 
various  filters  designed  at  8.00  and  5.33  KHz.,  designed  after  the 
specification  of  Table  1. 

2.4  Digital  Filter  Design 

The  filters  that  were  designed  for  the  MTI  application  under 
consideration  were  finite  impulse  response  filters,  i.e.  their  impulse 
response  was  non  zero  for  a finite  sample  number,  and  exactly  zero  for 
all  sample  times  greater  than  the  order  of  the  filter.  Filters  with 
finite  impulse  response  functions  exhibit  Gibbs  phenomena^  in  the 
frequency  domain,  i.e.  the  amplitude  vs.  frequency  response  have 
ripples  due  to  the  discontinuites  in  the  frequency  response  function, 
or  due  to  the  finite  digital  impulse  response  function^.  The 
filters  used  in  this  study  were  designed  with  the  Tchebeychef 
approximation  scheme  to  approximate  the  desired  pass  band  character- 
istics. Such  filters  have  the  property  that,  in  either  pass  or  stop 
bands,  amplitude  ripple  is  fixed  and  their  phase  is  linear.  Hence 
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the  terminology,  Tchebeychef  equiripple,  linear  phase,  finite  impulse 
response  filters, is  used.  Of  the  various  filter  designs  possible> 
the  equiripple  design  gives  the  minimum  transition  band  between  pass 
and  stop  band.  ^ illustration  of  the  frequency  response 

characteristics  of  such  filters  is  shown  in  Figure  1.  Table  1 
summarizes  the  properties  of  the  various  filters  that  were  designed 
and  tested  in  this  report.  Table  2 tabulates  the  actual  filters’ 
performance  as  a function  of  radial  velocity.  Table  3 lists  the  filter 
weights  used  in  this  report  whose  frequency  response  is  tabulated  in 
Table  2. 
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SECTION  III 


RAINFALL  CLUTTER  MODELLING  AND  MTI  OUTPUT 


3.0  MONTE  CARLO  SIMULATION  OF  CLUTTER  AND  TARGET 

The  performance  of  the  various  order  filters  was  evaluated  by 
simulating  the  amplitude,  phase,  and  velocity  of  both  rain  drops  and 
targets  by  means  of  a Monte  Carlo.  A simple  multiplication  of  the 
filter  power  spectrum  by  the  power  spectrum  of  rain  clutter 
return  signal  is  not  appropriate  in  unraveling  the  clutter's  re- 
jection in  the  presence  of  signal,  because  the  power  spectrum  of 
clutter  is  statistical  in  nature. 

The  complex  radar  return  from  an  assembly  of  M scatterers  and  a 
target  at  time  t is  given  by  Equation  8. 


e 


X 


(8) 


M * number  of  rain  drops  in  fiducial  volume 
<f>k  is  the  initial  phase  at  t = 0 of  the 


k th 

drop  (4^  = — — , R^  radial  distance  to  k drop) 

th. 

is  the  velocity  of  the  k drop. 


A^  is  the  amplitude  of  power  return  of  kth  drop. 


such  that 
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<1  E l\|2|  > = 20  dB 
k=l 


(9) 


o, j,  is  the  target  cross  section  and  is  such  that  <cr,j,>  = 0 dB> 
<()^  is  the  initial  phase  of  the  target  (4ttRj,/X),  where 
Rj,  is  the  radial  distance  of  the  target. 

Both  the  target's  initial  phase  and  the  initial  phase  of  the  scatter- 
ers  are  assumed  to  be  statistically  independent  and  uniformly  distri- 
buted, viz, 


p(<f>)  <4  = . 


(10) 


In  the  Rayleigh  limit  (D<<X)  0*-^ » l4]  the  ampiitude  must  vary 
the  cube  of  the  diameter.  The  diameter  follows  a Marshall-Palmer 
distribution 


as 


N (D)dD  = e"AD  d(DA)  (11a) 

A = 4.1  R-*21,  (lib) 

where  A is  in  units  of  mm  ^ and  R,  the  rain  rate,  is  in  units 
of  mm/hr. 
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3.1  Rain  Clutter  Cross  Section  Model 


The  (total  differential)  cross  section  of  rain  drops,  approximated 
by  dielectric  spheres, is  in  the  Rayleigh  limit  (diameter/wavelength  <<  1), 


where: 

a = radius  of  sphere 

X * wavelength  of  incident  radiation 

n = complex  index  of  refraction  of  dielectric 


e - dielectric  constant  ^ /l.333 
6 = conductivity  'V'  0 

Q 

w = 2irf  = 2it  — = angular  frequency  of  radiation. 


given  V ~)  by 


a(0)  4ira2  ~~  — | (1  + cos2  6)  • sin  9 (12) 


(14) 


Defining 


> 


(15) 
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it  is  found  that  |k|^  = .9275  at  20°C  at  a wavelength  of  3.21  cm, 

whereas  at  0°C,  |k|  2 = .9300  t14} 

An  average  diameter  for  heavy  rain  ^ ^ (50  mm  per  hour)  is 

(2) 

.1  cm  which  yields  an  average  total  back  scatter  cross  section  of 
—10  9 9 

3.45  x 10  nr.  In  a fiducial  volume  of  5 x 10  nr  there  are 

11  2 
3 x 10  rain  drops,  which  translates  to  a cross  section  of  100  m . 

(Fiducial  volume  1°  x 1°  x 1 us  at  10  KM) . Small  aircraft  cross  sections 

2 

considered  in  this  report  are  in  the  neighborhood  of  1 m . Thus  a 
clutter  to  signal  ratio  for  the  rain  and  the  aircraft  considered  is 
of  the  order  of  20  dB. 

3.2  Rain  Velocity  Distribution  Model 

The  velocity  spectrum  of  rain  drops  has  been  hypothesized  ^ to 
be  a Gaussianly  distributed  spectrum.  The  width  of  the  distribution 
is  ascribed  to  four  effects 

°v  = aShear  + °Beam  + aTurb.  + aFall’  (16) 


and  in  general 


I aShear  I < ^ °BeaiJ  > 1 0Turb . I > ^ °Fall  I 

Appendix  A presents  an  analysis  of  the  origin  as  an  order  of  magnitude 

m 

calculation  of  the  various  terms.  Suffice  it  to  say  that,  crv<2  in 
thunderstorms.  The  important  feature  of  (?gjiear  that  it  is  linearly 
proportional  to  the  range  of  observation  and  is,  for  near  horizontal 
scanning,  the  dominant  contribution  to  the  variance  of  the  Doppler  spectrum. 
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3.3  Filter  Operation 


Three  filters  of  order  n = 15,  19,  and  23  respectively  were  used 
as  weights  in  a Digitized  MTI  filter.  Since  in  the  system  design 
application  two  PRF's  are  to  be  employed  as  noted  in  Section  2.1, 
there  were  three  filters  implemented  at  each  of  the  PRFs 

Denoting  the  response  of  filter  i,  lii<6  as  F^,  the  response  of 
filter  i to  samples  of  the  composite  radar  return  signal  (I,  and  Q, 
or  Amplitude  and  Phase)  given  in  Equation  (8)  is 


The  index  £,,  in  referring  to  R(&),  denotes  the  sample  of  the  composite 

signal  at  time  AT  where  T is  the  sampling  time,  which  is  also 
s s 

the  PRI  (pulse  repetition  interval) . 

Thus , 


Denoting  W(k),  1 < k < n,  as  the  equi-spaced  sample  response 
of  the  digital  filter  in  the  frequency  domain,  then  the  finite  im- 
pulse response  function,  w(£),  used  as  weights  in  the  transversal 


(18) 


R(A)  = ^ e^T  ej 


(19) 
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filters  considered  in  this  report,  is  obtained  by  the  Discrete  Fourier 
Transform 


w(£)  - - E W(k)  e_j 

n k=l  n 


(20) 


The  inverse  transform  of  the  impulse  response  function  is 


W(k) 


T w(Z) 
0*1 


j 2iTik 
n 


(21) 


3.4  Interpretation  of  Filter  Response 

In  the  special  case  of  no  clutter,  = 0,  in  Equation  8 and  for  a 
target  cross  section  aT  = 1,  the  response  of  the  filter  F^  is 

n 

F,  = I £R0O  w.(A)|  (22a) 

1 Z= 1 


n 


. x J ,itVtT 
e-j^T  e^  ^ w^i) 


wi(£)  e 


• 4 v^T  A 
J T s 


(22b) 

(22c) 


Moreover,  this  is  the  amplitude  of  the  Digital  Filter's  response  at  a 
Doppler  frequency  2v/X,  as  Equation  22c  is  just  the  Discrete  Fourier 
Transform  of  the  filter's  weights. 
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SECTION  IV 

DESCRIPTION  OF  RESULTS 


4.0  RESULTS 
Program  Implementation 

The  filters  were  tested  with  two  sets  of  boundary  conditions. 

The  first  involved  generating  clutter  whose  Doppler  characteristics 
were  Gaussian  with  a mean  radial  velocity  of  10  (20  knots)  and 

whose  spectral  width  (standard  deviation)  was  1 m/sec.  The  clutter 
had  an  average  power  of  20  dB.  The  second  condition  had  the  same 
clutter  characteristics,  but  in  addition  had  a target  with  a power 
return  of  0 dB  and  a radial  velocity  of  50  m/sec  (100  knots),  cf. 
Equation  (8). 

Figures  3 and  4 illustrate  the  statistics  of  the  real  (I  channel) 
and  the  imaginary  (Q-Channel)  parts  of  the  sampled  received  signal. 

The  power  p returned  (sum  of  squares  of  I and  Q)  had  statistics 
that  are  exponentially  distributed 


df  _ - p/<p>  d£ 

dp  6 <p>  (23) 


where  <p>  was  20  dB. 

If,  as  in  the  second  case,  there  was  a target  present,  then  the 
average  return  clutter  to  signal  power  was  20  dB.  In  taking  the  square 
root  of  the  power  spectrum,  the  Rayleigh  spectrum  of  amplitude  is 


observed. 
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df 

dA 


(24) 


0A  - A2/<A2>  dA  . r~ 

= 2A  « — - , A = /p  . 

<A  > 

Figures  3 through  6 are  samples  of  I,  Q,  power  and  amplitude  drawn 
from  500  independent  ensembles  of  rain  scatterers  and  targets  described 
by  Equation  (8).  Both  the  power  spectrum  and  the  amplitude  spectrum 
are  displayed  in  Figures  5 and  6. 

Figures  7 to  12  are  8 KHz  filters  with  15,  19,  and  23  elements  .' 
Figures  13  to  18  are  5.33  KHz  filters  with  15,  19,  and  23  elements 
respectively.  Figures  7 through  18  alternate  the  boundary  condition 
of  clutter  only  with  clutter  plus  target. 

These  figures  are  histograms  with  the  vertical  scale  relating  the 
number  of  occurrences  and  the  horizontal  scale f the  amplitude  output 
of  the  filter.  Note,  if  clutter  were  absent  and  the  signal  were  in 
the  passband,  then  (since  aT  = 1)  the  output  of  the  filter  is  1 + 6p 
(cf  Figure  1 and  Table  1) , the  particular  value  depending  solely  upon 
the  given  Doppler  frequency. 

Conclusions  can  be  drawn  from  these  figures  regarding  the  pro- 
bability of  detection  when  an  aircraft  is  present  by  selecting  a 
threshold  level.  However,  when  the  threshold  is  chosen  and  no  target 
is  present,  false  alarm  probabilities  can  only  be  grossly  estimated 
due  to  the  limited  statistics.  The  information  to  be  learned  from  the 
figures  is  in  their  internal  comparisons,  i.e.,  increasing  the  filter 

elements  sharpens  the  response  (more  processing  time  means  better 
frequency  selectivity) . 
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SECTION  V 


PROCESSING  LIMITATIONS 


5.0  SIGNAL  CORRELATION 

Nathinson  and  Reilly^^have  obtained  data  in  severe  rain  storms 
(%  40  mm/hr)  that  tend  to  support  the  notion  of  Gaussianly  distributed 
spectra.  Appendix  B derives  an  auto  correlation  function,  with  the 
assumption  that  individual  scatterers  have  velocity  distributed 
Gaussianly. 

In  the  continuous  case  the  lagged  products  are  defined  as 


P 


II 


(t) 


fl(t)  I(t+T)  dt 
Jl(x)  I (t)  dt 


(25) 


(T)  . dt  _ 

^ Jl(t)  I(t)  dt 

(26) 

where 

I(t)  = Rg  (R(t) ) , 

(27a) 

and 

Q(t)  = Im  (R(t)) 

(27b) 

In  the  Discrete 

(sampled)  case,  the  lagged  products  are 

- 

defined  (with  k and  L multiples  of  the  PRI)  as: 


N 

£l(k)  I (k+L) 
k=l 


£l(k)  I(k) 
k=l 
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(29) 


N 

g lW  Q (k+L) 

%<T)  ' fXft)  w 


The  function  p(x)  is  also  defined  as 


p(x) 


/TTT 

P„(T) 


P?Q(T) 


It  is  proved  in  the  appendix  that  for  Gaussianly 
velocities  that: 


4ttt  a 2 

Pw  = .-i(-ri  . 


Pll<*>  - Cop  P (x)  , 

PIQ(x)  - Sin  PCX) 


(30) 


distributed 


(31) 


(32) 


(33) 
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Thus,  p (t)  is  the  correlation  function  of  the  fluctuations,  with 
the  average  Doppler  term  removed.  The  Fourier  transform  of  P(t) 
displays  the  fluctuation  spectrum: 


S (f) 


'/* 


-J2irfT  p(t)dTi 


(34a) 


S(f) 


(34b) 


Thus,  the  breadth  of  the  distribution  is  governed  by  2 o^/X,  viz, 
its  width  varies  linearly  with  the  variance  of  the  Doppler  spectrum 
and  inversely  proportional  to  the  wavelength.  Typical  values  of 
are  1 to  2 m/sec,  and  with  X = 3 cm,  the  width  of  the  fluctuation 
spectrum  is  66  or  133  Hz  for  ctv  of  1 or  2 respectively. 

Either  equation  32  or  33  will  allow  a determination  of  the 

4ttVt 

average  Doppler  velocity  shift  due  to  the  cosine  or  sine  term  of  — . 

A 

Figures  19  and  20  display  both  Pii(t)  and  p(x),  which  was  ob- 
tained from  the  data  of  J.  Roberts  Roberts  had  digitized  signals 

of  Nathanson's  experiment  on  collecting  I and  Q signals  of  rain  storms. 
The  radar  had  a wavelength  X = 5.2  cm,  and  a PRF  = 1.25  KHz.  Thus, 
from  Pj-j  (t)  , the  average  velocity  is  approximately  <V>^3.25  m/sec  = 

2 ^5 

6.5  knots,  and  the  width  ov  = < (V  - <V>)  )^  is  approximately  1 m/sec. 

With  a knowledge  of  the  correlation  time  determined  by  p(x),  an 
upper  limit  is  placed  upon  coherent  processing  of  signals  determined 
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by  the  widths  of  Doppler  velocity  spectra.  Or  stated  in  the  frequency 
domain,  the  width  of  the  Doppler  spectra  is  governed  solely  by  the 
width  of  the  velocity  spectra*  and  this  width  placed  the  upper  limit 
upon  a filter's  rejection  capabilities,  i.e.  the  amount  of  power  that 
is  within  the  pass  band  is  determined  not  only  by  the  central  Doppler 
frequency,  but  also  by  the  width  of  the  Doppler  spectrum  (as  is  expected). 

Using  Equation  31,  an  estimate  of  the  upper  limit  of  the  pro- 
cessing time  can  be  made  as  follows.  If  the  correlation  of  the  signal 
from  the  beginning  to  the  end  of  the  processing  is  .1,  then  for  a 
wavelength  of  3 cm,  and  a v = x = ■ /-2fl.n  .1  ■ 5.2  m sec.  This 

corresponds  to  a 41  element  filter  at  8 KHz  and  27  elements  at  5.33  KHz. 
5.1  Quantized  Noise 

The  process  of  sampling  and  digitizing  an  analog  signal  introduces 
noise  (random  error)  into  the  signal  that  is  being  analyzed.  The 
noise  is  due  to  the  quantization  of  the  analog  signal. 

Appendix  C discusses  the  necessary  conditions  for  the  noise  to  be 
considered  white  noise,  viz,  that  successive  signals  being  digitized  are 
uncorrelated  within  a quantum  level.  A derivation  of  the  ratio  of 
white  noise  due  to  quantization  of  the  true  signal  is  presented  in 
the  appendix.  The  ratio  of  quantized  noise  power  to  received  power 
for  a 10  bit  converter  is  -52  dB,  an  11  bit  converter  is  -58  dB,  etc. 
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SECTION  VI 


CONCLUSIONS 


Digital  filters  with  increasing  elements,  but  having  approximately 
the  same  frequency  response  characteristics,  perform  better  on  suppress- 
ing clutter  than  filters  with  a smaller  number  of  elements.  It  is  also 
of  importance  to  note  that  the  processing  time  should  not  be  greater 
than  the  time  determined  by  the  autocorrelation  of  the  clutter  to  be 
suppressed.  If  the  processing  time  exceeds  the  coherence  time  of 
clutter,  then  the  digital  filter  will  be  processing  parts  of  the  data 
stream  that  are  statistically  independent  of  the  preceding  portion; 
hence  the  filter  will  be  used  as  an  incoherent  device.  Stated  in 
context  of  an  X-Band  radar,  the  signal  processing  must  be  done  in  a 
time  of  the  order  of  5.0  milliseconds  if  the  spectral  width  of  the 
clutter  is  1 m/sec. 
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TABLE  1 


FILTER  DESIGN  CHARACTERISTICS 
All  Entries  Refer  to  Figure  1 


N 

fl 

f2 

f2-fi 

6P 

«s 

fl 

f2 

f2-fi 

6P 

«s 

15 

.398 

1.333 

.935 

1.03 

-48.02 

.670 

1.333 

.664 

1.04 

-47.94 

19 

.626 

1.333 

.707 

1.03 

-48.00 

.821 

1.333 

.502 

1.04 

-47.93 

23 

.7674 

1.333 

.576 

1.04 

-47.98 

.932 

1.333 

.401 

1.05 

-47.89 

8.0  KHz  Filters  5.33  KHz  Filters 


frequency  in  Kilohertz 

6_ , 6 in  dB 
p*  s 

N is  the  order  of  the  filter 
Note  V = Xf/2,  thus  1.333  KHz 

— @ X = 3 cm. 


corresponds  to  a Doppler  velocity  of  20 


sec 


TABLE  2 


Frequency  Reponse  of  Filters  (dB)  versus  Radial  Velocity  for  X = 3 cm 


8.00  KHz  5.33  KHz 


„ ORDER 

15 

19 

23 

15 

19 

23 

VR 

0 m/sec 

-48. 

-48. 

-48. 

-48 

-48. 

-54. 

10 

-17.7 

-37.1 

-48. 

-50.5 

-60 

-59 

20 

-.5 

-.5 

-.6 

-.6 

-.6 

-.6 

30 

-.4 

-.5 

+ .02 

-.03 

+ .5 

+25 

40 

+ .46 

+.02 

-.3 

-.5 

+ .5 

-.56 

50 

-.5 

+ .4 

-.5 

-.3 

+.46 

+ .25 

60 

+ .5 

-.5 

+ .02 

-.6 

-.6 

-.6 

70 

-.5 

+ .4 

-.4 

-50.5 

-60. 

-59 

80 

+.5 

+ .02 

-.2 

-48. 

-48. 

-48 

90 

-.3 

-.5 

+.02 

-50.5 

-60. 

-54 

100 

-.5 

-.5 

-.5 

-.6 

-.6 

-.6 

110 

-17.7 

-37.1 

-48 

-.03 

+ .5 

+.25 

120 

-48. 

-48. 

-48 

-.56 

+.5 

-.56 

130 

-17.7 

-37.1 

-48 

-.3 

+.46 

+.25 

140 

-.5 

-.5 

-.6 

-.6 

-.6 

-.6 

150 

-.4 

-.5 

+.02 

-50.5 

-60. 

-54. 

160 

+ .46 

+.02 

-.3 

-48. 

-48. 

-48. 

170 

-.5 

+0.4 

-.5 

-.6 

-.6 

-.6 

180 

.5 

-.5 

-.5 

-.6 

-.6 

-.6 

190 

-.5 

+.4 

-.4 

-.3 

+ .4 

+.  3 

200 

+ .5 

+ .02 

-.2 

-.6 

+ .5 

-.6 

210 

-.4 

-.5 

+ .02 

-.3 

+ .46 

+ . 25 

220 

-.5 

-.5 

-.5 

-.6 

-.6 

-.6 

230 

-17.7 

-37.1 

-48. 

-50.5 

-60. 

-54. 

240 

-48. 

-48. 

-48. 

-48. 

-48. 

-48. 

X = 3 cm.,f  is  Doppler  frequency  from  target  with  Radial  Velocity 

VD  . / 

R in  m/sec. 
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TABLE  3 


FILTER  WEIGHTS 


8.00  KHz  5.33  KHz 


n=15 

n=15 

W(1) 

— 005067 

.010351 

W(2) 

+.041367 

-.035465 

W(3) 

+.025585 

-.0044975 

W(4) 

-.00634 

+.05458 

W(5) 

-.069283 

+.05698 

W(6) 

-.14857 

-.082521 

W(7) 

-.21513 

-.2988 

W(8) 

+.75885 

.59474 

Where  terms  greater  than  8 have  been 
by  W(16-k)  = W(k),  I s k i 8. 

omitted  and  are 
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TABLE  3 (Cont'd) 


8.00  KHz; 

5.33  KHz 

n=19 

n=19 

W(1) 

-.0097096 

-.011098 

W(2) 

-.027526 

.029622 

W(3) 

.021661 

.000133 

W(4) 

.033945 

-.035468 

W(5) 

.043443 

-.021503 

W(6) 

.012459 

.0524 

W(7) 

-.05868 

.075366 

W(8) 

-.15262 

-.066455 

W(9) 

-.23279 

-.30694 

W(10) 

+.73566 

.57191 

W(20-k)  = W(k),  1 < k < 10 
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TABLE  3 (Cont’d) 


8.00  KHz 

5.33  KHz 

n«23 

n=23 

W(1) 

.023511 

.011699 

W(2) 

-.016915 

-.026468 

W(3) 

-.022938 

.0013783 

W(4) 

-.017912 

.025923 

W(5) 

.0049624 

.010138 

W(6) 

.036139 

-.036618 

W(7) 

.051167 

-.033569 

W(8) 

.024877 

.047031 

W(9) 

-.050044 

.085507 

W(10) 

-.15301 

-.054856 

W(ll) 

-.24261 

-.31104 

W(12) 

.7219 

+.55771 

W(24-k)  = W(k) , 1 < k < 12 
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|i*-43,2s7] 


f2-f,  = TRANSITION  WIDTH 

Ts  = SAMPLING  TIME 

l./Ts  s SAMPLING  FREQUENCY 

Ss  = MAXIMUM  DEVIATION  IN  STOP  BAND 

Sp  = MAXIMUM  DEVIATION  IN  PASS  BAND 

f,  = Ts--  f2 


Figure  I CHARACTERISTIC 


OF  EQUIRIPPLE  DIGITAL  FILTER 


AMPLITUDE  (dB)  $ 0 RESPONSE  (dB) 


5.33  KHz  I0.66KHZ  I6.0KHZ 


VELOCITY  (m/sec) 


Figure  2 AMPLITUDE  RESPONSES  OF  8.00  KHz  AND  5.33  KHz  FILTERS 
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FIGURE  3 STATISTICS  OF  I CHANNEL 
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FIGURE  4 STATISTICS  OF  Q CHANNEL 
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FIGURE  5 STATISTICS  OF  RECEIVED  POWER 
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FIGURE  6 STATISTICS  OF  RECEIVED  AMPLITUDE 
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FIGURE  7 OUTPUT  STATISTICS  OF  15  ELEMENT  FILTER  WITH  PRF  OF  8Khz,  CLUTTER  ONLY 
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FIGURE  8 OUTPUT  STATISTICS  OF  15  ELEMENT  FILTER  WITH  PRF  OF  8Khz,  TARGET  PLUS  CLUTTER 
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FIGURE  9 OUTPUT  STATISTICS  OF  19  ELEMENT  FILTER  WITH  PRF  OF  8.0  Khz,  CLUTTER  ONLY 
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FIGURE  10  OUTPUT  STATISTICS  OF  19  ELEMENT  FILTER  WITH  PRF  OF  8.0Khz,  TARGET  PLUS  CLUTTER 
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FIGURE  11  OUTPUT  STATISTICS  OF  23  ELEMENT  FILTER  WITH  PRF  OF  8.0Khz,  CLUTTER  ONLY 
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FIGURE  12  OUTPUT  STATISTICS  OF  23  ELEMENT  FILTER  WITH  PRF  OF  8.0Khz,  TARGET  PLUS  CLUTTER 
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FIGURE  13  OUTPUT  STATISTICS  OF  15  ELEMENT  FILTER  WITH  PRF  OF  5.33  Khz,  CLUTTER  ONLY 
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FIGURE  14  OUTPUT  STATISTICS  OF  15  ELEMENT  FILTER  WITH  PRF  OF  5.33  Khz,  TARGET  PLUS  CLUTTER 
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FIGURE  15  OUTPUT  STATISTICS  OF  19  ELEMENT  FILTER  WITH  PRF  OF  5.33  Khz,  CLUTTER  ONLY 
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FIGURE  16  OUTPUT  STATISTICS  OF  19  ELEMENT  FILTER  WITH  PRF  OF  5.33  Khz,  TARGET  PLUS  CLUTTER 
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FIGURE  17  OUTPUT  STATISTICS  OF  23  ELEMENT  FILTER  WITH  PRF  OF  5.33Khz,  CLUTTER  ONLY 
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FIGURE  18  OUTPUT  STATISTICS  OF  23  ELEMENT  FILTER  WITH  PRF  OF  5.33Khz,  TARGET  PLUS  CLUTTER 
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Figure  19  RAIN  CLUTTER  AUTOCORELLATION  puU) 
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p( t)  correlation  function 


Figure  20  RAIN  CLUTTER  CORRELATION  p{ r) 


APPENDIX  A 

RAIN  CLUTTER  SPECTRAL  DISTRIBUTION 


Consider  a radar  that  has  an  illumination  function  defined  as 
l(0,ij>),  such  that 


Also,  assume  that  <}>  and  0 are  independent  and  Gaussianly  distributed. 


If  a radial  velocity  can  now  be  formulated  in  terms  of  either 
or  both  of  these  angles,  then  the  spectrum  due  to  wind  shear  and 
beam  broadening  can  be  computed. 

For  wind  shear,  the  velocity  as  a function  of  elevation  height 

h is 


1(6, 4)  d0d<j>  - 1 


(1) 


(2) 


V(h)  - VQ  + k h , 


(3) 


W = k 


(4) 
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In  the  event  that  the  elevation  angle  is  small,  then  the  height 


h can  be  related  to  the  range  and  elevation  angle  as 

h = R sinG  % R6  , (5) 

= V + kR9  (6) 

o 

Assuming  that  the  intensity  pattern  in  the  elevation  angle  is 
distributed  as 

(9-9o)2 

- -I  . 

, 2a  * 

p(e)d9  = r (7) 

/2ir  Og 

then  the  average  velocity  and  the  standard  deviation  are 


<v>  = 

p(0)  v(0)d0 

(8) 

<v>  = 

V + k R 0 
o o 

(9) 

Q 

< to 
II 

<(V-<v>)^  = R^ 

(10) 

0 = 
V 

k R o0 

(ID 
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The  important  feature  of  this  result  is  that  wind  shear  causes 
the  Doppler  velocity's  spectrum  width  to  increase  linearly  with  range. 

Beam  broadening  is  the  term  that  applies  to  wind  driven  clutter 
perpendicular  to  the  bore  sight. 

Since  the  illumination  pattern  subtends  other  sectors  than  the 
boresight  the  radial  velocity  is 


where  V is  the  velocity  of  the  wind  perpendicular  to  the  boresight 
In  the  case  of  narrow  gaussian  beams 


or  Beam  Broadening  is  dependent  linearly  on  the  cross  wind  velocity. 

Turbulence  describes  the  statistical  properties  of  the  wind, 
and  as  such  describes  the  variation  of  wind  speed  independent  of 
height,  range,  etc.  As  such,  it  has  been  found  that 


VR(+)  - V sin  ♦ 


(12) 


"V  = 0 


(13) 


(14) 


TURB 


[17,18] 


(15) 


a 


v 
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Fall  velocity  contributes  to  the  spectrum  negligibly 


vert  lm 

a 

v sec 


(16) 


Moreover  for  low  elevation  angles 


a = a Sin0  'v*  0<j  « — 

v v v sec 


(17) 


which  is  entirely  negligible. 

Beam  broadening  for  = 1° corresponds  to  'v*  V cross  wind. 
If  the  cross  wind  velocity  is  20  — — , then  for  Beam  Broadening 

S 6 C 


«4tn 

av  sec 


(18) 


Shear  constant  k = W has  been  found  ,*-4  to  be  approximately  5.7 
thus  at  a range  of  10  km,  and  for  a 1°  Beam 


m 

Sec*Km 


Shear  - _ 5.7  lm 

a = 10  x -r=~  ~ 

v 57  sec 


(19) 


It  is  not  implausible  to  believe  that  these  distribution  widths 
are  independent,  hence 


2 „ a2TUKB+  c2FALL  + a2B.B  + a2SHEAR 


v V 


(20) 
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APPENDIX  B 


FLUCTUATION  CORRELATION 

The  autocorrelation  functions  denoted  byp^j(i),  Pj-q(t)  and  p(t) 
are  experimentally  averaged  with  respect  to  time. 

There  are  theorems  [20] 

proved  to  state  the  ensemble  averages 
equal  averages  with  respect  to  time  for  stationary  processes . 


+T 

T^°  Tir  -T  S(t+t)dt  = < f(t)  g(t+x)  > , 


(Bl) 


where  the  brackets  denote  a statistical  average  over  the  ensemble, 
and  the  time  t is  arbitrary, as  are  the  function  f(t)  and  g(t),  with 
the  provision  that 


<f(t1)  = <f(t2)>  (B2) 

<g(t1)  = <g(t2)>  (B3) 

as  well  as  all  the  moments  of  f and  g. 

In  the  present  case,  the  function 

N 4 V t 

f(t)  = E \ cos  (<f>k  + -yL-)»  - (B4> 

k=l 
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(B5) 


N 

g(t)  - £ \ cos  + 
k=l 


4irVkt 

T~ 


) 


or 


v.  4uV,t 

•^2  sin  (<j>k  + — ^ — ) (B6) 

k=l 


Define 


, , A 5Re  R(t)  Re  R(t  + t)  dt 
PII(T)  = J Re  R(+)  Re  R(t)  dt 


(B7) 


where  Rg  denotes  real  part. 


A JR£  R(t)  Im  R(t+x)  dt 
PIQ(T)  = J Rg  R(t)  Re  R(t)  dt 


(B8) 


where 


R(t)  = 


N 

L *k 

K=1 


e^  (<j>. 


4ttV,  t 

k + -^> 


(B9) 


describes  the  and  channel  of  the  received  signal  from 
N scatterers  at  time  t,  each  of  which  has  an  initial  phase  <J>,  , and 

4TrVk 

doppler  velocity  — — 
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Lemma  I 


For  Gaussian  velocity  distributions 

(v-v)2 

2o2 

h(v)dv  = — 

m ov 


(BIO) 


j*Re  R(t)  R(t+x)dt  = !<A2>  e 


0_2  2 2 
8tt  t a 

v 


X ~ X2 


(B11) 


. , v 4itVt 

then  Pjj(t)  = cos  ~ — 


o 2 2 2 

8it  t 

V 

, 2 


(B12) 


, , „ , 4ttVt 

and  Pjq  (t)  = sin  — — 


8it2t2o  2 
v 


(B13) 


Thus 


P (t)  = e 


o 2 2 2 

8n  T o 

V 


(B14) 


Therefore, 


R(t)  R(t  + t)  dt 


<R  R(t)R(t  + t)> 
e 


(B15) 


X 


(t+T 


> 


(B16) 
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where  the  average  is  performed  over  A^>  4^  and  . 

Assuming  that  the  drops  are  statistically  independent, 
use  of  the  following  identities 


Cos  A Cos  B = y Jcos  (A  + B)  + Cos  (A  - B)J 

» 

Cos  A Sin  B = j |sin  (A  + B)  + Sin  (A  - B)| 


and 


<cos  <j>,  > = < sin  <|>,  > =0,  then 


with 


Airy,,  t 
k T X 


471 V,, 

A = <k.  + - k 


B = 6 + ilV  <t+T) 

X X 


<cos  (*k  + iiat  )aJ  «(  + 4^  (t+T)>  . 1 ej  ( 4^  } 4 


2 *k  X > k£ 


whence  the  average  over  velocity  and  amplitude  give 


fN  2 \ 4 V x 

< Re  R(t)  R(t  + r)dt  = j <A  > < eJ  > 


then  making 

(B17) 

(B18) 

(B19) 

(B20) 

(B21) 

(B22) 

(B23) 
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Finally,  if  Gaussian  velocity  distribution  is  assumed 


<Re  R(t)  R(t+x) 


N 1 4nVt  - 8»2t2c  2 

2 <A  ” e T e — v 


Therefore 


pn  (T> 


Cos 


4ttVt 

X 


-8tt 

e 


2 2 
T 


Sin 


4ttVt 

X 


2 2 
-8tt  x 


and 


P (x) 


0 2 2 
- 8tt  t 

e 7 


The  power  spectral  density  is  also  easily  found  to  be 


but 


3(f)  = Jej2lTfTjR  (t)  R*  (t  + t)  dT  , 
S(f)  = < R(t)  R*  (t+  x)  > dx  , 


(B24) 


(B25) 


(B26) 


(B27) 


(B29) 
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with 


*00  »*(«h)  - < [E  (\  e%  + .-Jfl  *4fl 


N <A2>  < eJ  ‘"-p  > 


„ ..2.  . -1  4r  V T -8n2  T2  cr  2 

N <A  > > eJ  — - — e — v 

A X 


where  a Gaussian  distribution  for  velocity  is  assumed. 
Thus 


S(f) 


j*eJ2lrfT<  R(t)  R*  (t+x)  > dx 

Jj2irfx  „ .2^  -j4irV  x -8tt2  x a 2 
eJ  N <A  > eJ  y e — j v 


Noting  that 


per  - t*/2«2  = 


v^T 


with 


a 


and  8 = 2irf  - 


4ttV 

X 


* 


then 


S(f)  * 


N <A2>  X 
4tt  a 

v 


- x2 

e 8ir  2 
v 


(f 


2V.2 

X; 


,<H]> 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 

(B35) 

(B36) 

(B37) 
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APPENDIX  c 
Quantization  Noise 

In  the  process  of  performing  analog  to  digital  conversion  of 

the  I*"*1  and  channel,  noise  is  introduced  into  the  filtering 

process  affecting  the  filter's  performance. 

In  sampling  the  1^  and  Q1”*1  channel,  if  the  time  between 

successive  samples  is  so  rapid  that,  for  a quantization  level  g,  the 

th 

autocorrelation  function  p (t) , (same  for  Q channel) , cannot  change 
appreciably  in  level,  then  systematic  biases  are  introduced,  i.e., 


1 “ Pjj(t)  g/2. 


(Cl) 


We  shall  not  consider  this  case,  but  rather  the  converse. 

The  text  showed  that  half  of  the  power  is  shared  by  the  I 
.th 


th 


and  Q channel  respectively.  If  each  channel  has  M bits  (1  bit 
for  sign),  and  if  the  digitizer  is  to  span  vo  deviation  of  the  voltage 
spectrum  of  I and  Q,  then 


_m-l 

g • 2 = vo 


(C2) 


Straightforward  probability  arguments  prove  that  if  the  signal 
populates  the  levels  uniformly,  with  G,  the  true  voltage  in  the  quantum 
interval,  then 


<G>  = 0 


(C33 


a2  = < (G  - <G>)  2> 


(C4) 


63 


(C5) 


(C6) 


12  • 4 


In  the  case  where  clutter  dominates  the  power  spectrum 


o2  = 1/2  £ |Aj2 


(C7) 


(the  factor  1/2,  as  the  power  divides  between  I and  Q) 

If  1 - Pjj(t)  > g/2  then  the  sampling  will  be 
independent,  and  the  noise  spectrum  white. 

The  fraction  of  white  noise  due  to  sampling  passed  by  the 
digital  filter  is  now 


2 f <G2>  |F(w)|2  dw 
5 dw 


(C8) 


For  the  digital  filters  in  this  study 


f 


@ 8.0  KHz 


@ 5.33  KHz. 


(C9) 
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Thus,  an  upper  limit  can  be  placed  upon  the  amount  of  noise  due  to 
quantization  passed  through  the  filters;  it  is 


2 <G2> 


2 

> <G  > 

'Xj 


i ,2  2 

1 va 

12  4m-l 


(CIO) 


(Cll) 


For  clutter  with 


°2=^l\|2 


(C12) 


= 50, 


v = 3 standard  deviation,  and  a 6 bit  converter,  the  power  passed  is 


32 .50 

12. 4b 


% .075 


(C13) 


th 


Thus,  the  ratio  of  quantized  noise  to  power  in  either  the  I or  Q 
channel  is 


th 


.075 

50 


= .0015,  or  28  dB 


(C14) 


For  an  8 bit  converter, the  ratio  is  40.0  dB. 
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